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Abstract 

The stochastic driving force exerted by a single molecular motor (e.g., a 
kinesin, or myosin) moving on a periodic molecular track (microtubule, actin 
filament, etc.) is discussed from a general viewpoint open to experimental test. 
An elementary "barometric" relation for the driving force is introduced that (i) 
applies to a range of kinetic and stochastic models, (ii) is consistent with more 
elaborate expressions entailing explicit representations of externally applied 
loads and, (iii) sufficiently close to thermal equilibrium, satisfies an Einstein- 
type relation in terms of the velocity and diffusion coefficient of the (load-free) 
motor. Even in the simplest two-state models, the velocity-vs.-load plots ex- 
hibit a variety of contrasting shapes (including nonmonotonic behavior). Pre- 
viously suggested bounds on the driving force are shown to be inapplicable in 
general by analyzing discrete jump models with waiting time distributions. 



Molecular motors are protein molecules such as myosin, kinesin, dynein, and RNA 
polymerase, that move along linear tracks (actin filaments, microtubules, DNA) and 
perform tasks vital to the life of the organism — muscle contraction, cell division, 
intracellular transport, and genomic transcription (1-5). Understanding how they 
operate represents a significant challenge. The hydrolysis of adenosine triphosphate 
(ATP), with the release of adenosine diphosphate (ADP) and inorganic phosphate 
(Pi), is known to be the power source for many motor proteins. An activated motor 
may well be in a dynamical or, better, a stochastic steady state but it cannot be in 
full thermal equilibrium. 

Striking in vitro experiments observing individual motor proteins moving under 
controlled external loads (6-11) have stimulated enhanced theoretical work aimed 
at understanding the mechanisms by which a biological motor functions. From a 
broad theoretical perspective, a molecular motor is a microscopic object that moves 
predominantly in one direction along a "polarized" one-dimensional periodic struc- 
ture, namely, the molecular track (1-11). In recent years, in addition to traditional 
chemical kinetic descriptions (see, e.g., (12) and references therein) and various more 
detailed schemes (11,13,14), so-called "thermal ratchet" models have been proposed 
to account for the mechanics: sec the review (15). 

A common feature of most approaches is that a motor protein molecule is asso- 
ciated with a labeled site I {— 0, ±1, ±2, . . .) on the track and is pictured as being 
in one of N essentially discrete states j, which may be free of (say, j = 0), or bound 
to ATP and its various hydrolysis products {j = 1, 2, . . . , — 1). Thus, for a ki- 
nesin molecule, K, on a microtubule, M, the (A^ = 4) states identified might be M-K, 
M-K-ATP, M-K-ADP-Pi, and M-K-ADP (8,12). Transition rates between these states 
can be introduced via 

Ui U2 Un-1 Um 

0^ # 1^ # ••• # (7V-1), - 0,+i, (1) 

Wi W2 Wn-1 Wn 

where the subscripts indicate that the states j are associated with successive sites, 
I and I + 1, on the track spaced at distances Ax — xj+i — xi — d : this defines the 
step size d. Of course, states ji, . . . , differ physically only in their spatial 
displacements d, 2d, . . . , nd, along the track. By the same token, the rates Uj and wj 
are independent of / (or x = Id); however, in the subsequent developments it proves 
useful to allow for spatially dependent rates Uj{l) and 'Wj{l). 

To properly represent physicochemical reality (that is, microscopic reversibility) 
none of the forward rates, Uj, or backward rates Wj may strictly vanish even though 
some, such as the last reverse rate, wn, might be extremely small (11,12). On the 
other hand, if, as one observes in the presence of free ATP, the motor moves under 
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no external load to the right (increasing x), the transition rates cannot (all) satisfy 
the usual conditions of detailed balance that would characterize thermal equilibrium 
if Eq. 1 were regarded as a set of chemical reactions (near equilibrium) between 
effective species ji (15). (Notice that one may envisage a second-order rate process, 
e.g., M-K + ATP ^ M-K-ATP, to conclude ui = A;i[ATP]; this can then lead to 
Michaelis-Menten type rate-vs.-concentration relations (6). However, one might also 
contemplate a small "spontaneous" or first-order background rate, Mi o > 0, that 
exists even in the absence of ATP.) 

Now, within statistical physics, the kinetic scheme in 
Eq. 1 represents a one-dimensional hopping process of a particle on a periodic but, in 
general, asymmetric lattice. After initial transients, the particle will move (16) with 
steady (mean) velocity V, and diffuse (with respect to the mean position, x — Vt, at 
time t) with a diffusion constant D. Complicated, but exact equations for V and D in 
terms of uj and wj have been obtained for all N (16), as exhibited in the Appendix. 
A dimensionless, overall rate factor which, rather naturally, appears (see Eq. Al), is 
given by the product 

This will play an important role in our discussion. Note, indeed, that viewing Eq. 1 
as a standard set of chemical reactions and requiring detailed balance would impose 
r = 1 (or £ = 0) whereas F > 1 (or e > 0) is needed for a positive velocity V. 
(One might comment, however, (17) that as regards the full chemistry, the complex 
of motor protein plus track may be regarded simply as catalyzing the hydrolysis 
of ATP: the reaction rates for this overall process may then be expected to satisfy 
detailed balance.) 

The simplest or "minimal" physical models have N = 2, and one can then cal- 
culate analytically not only the steady state behavior but also the full transient 
responses, specifically, the probabilities, Pj{l;t), of being in state ji ("at" site /) at 
time t having started, say, at site Z = in state j — time t — 0: see (17). In 
Ref. (17) only the special (limiting) cases with = W2 — were treated; but as 
seen below, this limit can be misleading and so for completeness (and for possible 
comparisons with experiment and simulation), we present the general N = 2 results: 
see Appendix. In particular, the velocity and diffusion constant for N = 2 are 

V = (^^^^-^^^^^d ^ (r - l)u;d, (3) 

Ui + U2 + Wi + W2 

D = 1[t + 1-2{T -lfuj/a]ujd\ (4) 

where P = U1U2/W1W2 (as in Eq. 2) and, for convenience, we have introduced the 
associated overall rates 

a = ui + U2 + wi + 'W2, uj = wiW2la. (5) 
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One can, of course, envisage more complicated schemes than Eq. 1, with various 
internal loops, parallel pathways, etcQ In all cases, however, there will be a well 
defined (zero-load) steady-state velocity V and a diff'usion constant D (independent 
of the particular states, j) which are susceptible to estimation by simulation even 
though their explicit mathematical expressions may be intractable. Furthermore, in 
real systems both V, as often demonstrated (6,8-11), and D (7) are susceptible to 
experimental measurement. 

Now there arises an obvious but crucial question, namely: "What (mean) driving 
force, /, will such a general motor protein model exert as it moves along the track?" 
That is the issue we address here from a theoretical standpoint. 

ANALYSIS 

Maximum Driving Force. The hydrolysis of one ATP molecule releases a 
free energy 

AGo of about 0.50 x 10~^^J (corresponding to 7.3 Kcal/M or 12 ksT at 
typical in vitro temperatures, T (3)). If all this free energy could be converted into 
mechanical energy and move the motor protein through a distance Ax = d, the step 
size (going from state 0^ to O^+i: see Eq. 1), the force exerted would be 

U,, = AGo/d. (6) 

Accepting that one molecule of ATP is sufficient to translocate the motor protein by 
one step (7), this expression clearly represents the maximal driving force that can be 
exerted. For a kinesin moving on a microtubule (6-12) with d ~ 8.2 nm (9) it yields 
fmax — 6.2 pN. If / is the driving force actually realized, the efficiency of a motor 
protein may sensibly be defined by S = f / fmax- 

Einstein Force Scale. Consider a small ( "mesoscopic" ) particle with "instan- 
taneous" position x{t) and velocity v{t) that undergoes one-dimensional Brownian 
motion in a fixed, slowly varying external potential, $(x). Under a constant external 
force, F = —{d^/dx), the particle diffuses with a diffusion constant which, for long 
times, t, satisfies 

/2t, (7) 

where ( ■ ) denotes a statistical average (18-21). In addition, the particle experiences 
an (effective) frictional force, fs = C'^{t), where C is a friction coefficient determined 
by the environment (18-21). In a steady state, the friction balances the external force, 

* Thus a backwards reaction directly from, say, state jj to 0; could account for "futile" ATP 
hydrolysis, i.e., without forward motion of the motor (13); but note that within N ^ 2 models 
(which enforce = 1) the phenomenon may be described simply by including the futile-hydrolysis 
parallel reaction rate in the backward rate wi . 
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F, leading to a drift motion, {x{t)) ~ Vt, with mean velocity V = F/( = Je/C- Now, 
by definition, Brownian motion takes place within full thermal equilibrium: that fact 
dictates (18-21) the Einstein relation ( — ksT/D which, in turn, implies the result 

fE = ksTV/D. (8) 

In the present context this is an appealing formula in that it sets a force scale in 
terms only of the velocity, V ^ and the diffusion constant, D, predicted by a motor- 
protein model (or observed in an experiment or simulation); one might call it the 
"Einstein scale" . However, because an activated molecular motor is not a Brownian 
particle and cannot be described by thermal equilibrium, there are no grounds for 
expecting to be related to the proper driving force, /. Nevertheless, we will 
show that in a certain limit such a Brownian motion "mimic" of an activated motor 
protein does provide an appropriate prediction for /. Indeed, Ref. 17 accepted the 
identification / = /e without discussion and used Eq. 8 to estimate driving forces for 
restricted {w2 — 0) N — 2 models: the values of / so obtained were not unreasonable 
in comparison with experiments (17). 

It should, perhaps, be mentioned in passing that Ref. (22) (see also (13)) invokes 
an Einstein relation in an analysis of "protein friction". However, this is a rather 
different context in which many "blocked" motor proteins are attached to a substrate 
and a rigid microtubule diffuses, apparently freely, close-by in the medium above. 
Quantitative arguments (13,22) explain the large frictional slowdown observed (rela- 
tive to an Einstein-relation estimate) as due to weak protein binding and unbinding. 

Barometric Formulation. Although the identification of the motor driving 
force / with /g is unjustified, it would be desirable to have a soundly based, general 
expression for / which, like /e, does not entail any intrinsic modifications or exten- 
sions of the motor model or of the associated physicochemical picture. To that end, 
consider the placement of an "impassable block" or barrier on the molecular track, 
say, between sites L and L + 1 (^ 1) or at distance x = D = l^d from the origin 
x = (fixed, as wc suppose, by where the motor starts). Such a barrier may be 
realized theoretically by decreeing that all states ji ioi I > L -\- 1 are inaccessible; 
this may be achieved simply by setting one of the local forward rate constants, say, 
— L), equal to zero. No other rate constants need be modified; but if, perhaps 
to take cognizance of some aspects of a realistic barrier, further nearby rate constants 
are changed, it will have no consequences for the main conclusions. 

It is intuitively clear that running a molecular motor up to such a barrier will 
lead (provided it does not detach from the track or "freeze" irreversibly, as might 
happen in practice (6,10)) to a stationary probability distribution, Pj{l,t ^oo) = 
P^{L — I), in which z = {L — l)d = D — x measures the distance back from the 
barrier. On very general theoretical grounds one should expect this distribution to 
decay exponentially with increasing z (except for possible deviations close to the 
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barrier) so that, explicitly, one has 

P°°{z/d) ^ Aje-"'. (9) 

The (positive) decay constant k, should, in principle, be experimentally measurable 
(although this may be difficult if nd is large). The amphtude ratios Aj/Aq must 
depend on the various rate ratios, Ui/wi, while Aq is set simply by normalization. 

To justify this surmise for the kinetic scheme in Eq. 1 (although it is of general 
validity), note that the mean flow between adjacent states 0/ and {N — l)/_i and 
between ji and (j — 1); [for j = 1,2, . . . , {N — 1)] must vanish for a stationary 
distribution. Balancing local forward and backward rates thus yields 



u,{l)P°Z,{L-l) = w,{l)P^{L-l), (10) 

for j — 1,2, . . . ,{N — 1). Starting from an initial nonzero value Pj°(0), one can 
then recursively determine P,^i(0), P^2(0), ■ ■ ■ , ^o°°(0), P^_i(l), ^^-2(1), ' ' ' • By 
induction, this leads directly to Eq. 9 [since the Uj{l), and 'Wj{l) become independent 
of / for, say, I < L — Iq where Iq is a fixed integer]. Most crucially one finds (with the 
notation of Eq. 2) that the decay constant is simply 

K^{\nT)/d^e/d. (11) 

Now, to interpret these results, consider a dilute gas of molecules of mass m 
moving in a gravitational field that acts "downwards" along the z axis. Each molecule 
then has a weight fa — fng; in addition, the equilibrium density distribution is given 
by (23) 

p{z) = p(0)e-"^^"/'^«^, (12) 

where p(0) is the density at the level z = 0. (Any deviations arising close to the 
"lower" wall (at 2; ~ 0) due to molecular size, structure, etc., have been neglected.) 
Comparing this well known barometric formula with the distribution Eq. 9 leads us 
to identify the driving force / of the molecular motor with 

fn = ksTiln T)/d = kBTe/d. (13) 

The subscript B here serves merely to indicate the barometric analogy underlying 
our identification. By comparison with Eq. 6 for fmax, we may expect e < AGq/^bP 
for a real molecular motor. 

Before studying this result in relation to extensions of Eq. 1 needed to describe 
a motor functioning under external loads, let us compare with /g. 
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Barometric vs. Einstein Scale. Suppose the molecular motor operates close 
to equilibrium in the sense that £ = In F is small. Then, on expanding in e at fixed 
a;/(7, Eqs. 2-5 and 13 yield 

JbUe = 1 + [i^ - {i^l<y)V - \{uj/a)s' + • • • , (14) 

for N = 2. Evidently, the coefficient of £ vanishes identically! Furthermore, one 
finds < u/a < so that the coefficient of e"^ is small, lying between ^ and ^. 
Consequently, and as anticipated, the Einstein scale approximates the barometric 
result very well when the motor operates sufficiently close to equilibrium. Indeed, for 
F < 10, calculations show that /b can exceed /e by no more than 44%. Furthermore, 
the series truncated at 0(£^) in Eq. 14 is reasonably accurate up to e ~ 5 (F ~ 150) 
where one has 1.473 < Jb/ Je < 2.535; beyond that, the bounds < fn/ Je ^ 
are effective. 

These specific results are hmited to N — 2; but we suspect (and have checked for 
N — 3) that the vanishing of the 0(e) term in Eq. 14 is independent of N. Likewise, 
we expect fs always to rise steadily above Je when e increases. Indeed, on recalling 
Eq. 2 for F, one observes from Eq. 13 that J'b is unbounded above and so, with an 
injudicious assignment of rate constants, it may even exceed f^ax (as given in Eq. 1)! 
Conversely, one may show from Eqs. 3, 4 and 8, that Je for N — 2 is bounded above 
by AksT/d (17). However, we will demonstrate below that this bound on /e is rather 
artificial and does not apply for models that account directly for the discreteness of 
ATP hydrolysis. 

Stalling Force Measured by Spring Compression. In a typical experiment 
on motor proteins (6-10), optical tweezers are used to carry a silica bead coated with 
a few molecules of the motor protein up to the molecular track. Then a single motor 
binds to the track and starts to move, exerting a force against the opposing load, F, 
as it pulls the bead towards a side of the optical trap. The external force F is a linear 
function of the displacement from the trap center, and the constant of proportionality 
can be measured. Thus the trap and bead work like a calibrated spring acting on the 
molecular motor. To represent such an experiment the load-free scheme embodied in 
Eq. 1 must, clearly, be extended. 

To this end, suppose the motor moves in a slowly varying external potential, ^{x), 
so that in translocating from site / to / + 1, additional mechanical work A$(a; = 
Id) — ^{x + d) —^{x) must be done (relative to tlic load- free situation). Of course, 
this corresponds to imposition of a local external force, F{x) — A^{x)/d, directed 
negatively. For an (ideal) optical trap of spring constant K we may take 

^x) = \Kx^, F{x) = K{x + Id). (15) 

In such a situation the motor should, in effect, compress the spring and, as t 
increases, attain a stationary distribution, say Pq{1), where, for simplicity, we focus 
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only on the (free) states 0;. This distribution should peak at some Is, corresponding 
to a mean (or most probable) compression of the spring by a displacement xs = Isd- 
Then the measured "stalling force" would be fs — Kxg- 

Now it is physically clear that under any local load, F(x), the transition rates, 
Uj{l) and Wj{l), must change. If, as traditional, one views the chemical transitions 
between successive states, j and j + 1, as proceeding in quasiequilibrium over various 
free energy barriers (13), one expects (in leading approximation) the rates to change 
exponentially with F{x)d/kBT. But how the exponential loading factors should be 
distributed among the various reaction processes, j ^ (j + 1), is far from clear: indeed, 
this distribution is of considerable interest in understanding the motor mechanism at 
a microscopic level. Without prejudice, therefore, we will explore the quasiequilibrium 
hypothesis that under a local load, F, the local transition rates change in accord with 

w.^ njf^ wfe+'i^'/''-^. (16) 

The distribution factors, Oj' and 9~ , need not be of uniform sign: but we certainly 
expect the overall factor 

N 

^-E(^;+^7)' (17) 

to be positive, implying an opposition to motion. Indeed, should the motor undergo 
diffusion in thermal equilibrium when not activated by ATP (as suggested parenthet- 
ically in the introductory discussion of Eq. 1), detailed balance considerations would 
dictate 9 — 1. As a supplement to our quasiequilibrium hypothesis this value of 9 
is also plausible for an activated motor that operates not too far from equilibrium. 
Notice that a negative 9^ or 9j simply means that the corresponding forward rate, 
uj, is enhanced, or the reverse rate, wj, is diminishedhy the internal molecular strain 
induced in the motor by the load. 

Accepting Eq. 16 we can find the stationary distribution Pq{1) with the aid of 
the rate-balance Eqs. 10 [ replacing P°°{L — I) by Pf{l) and the rates Uj and Wj 
in accord with Eq. 16 ]. The most probable motor location. Is, follows by equating 
P(f (/) and P(f (/ + 1) which leads directly to the condition 

N 

r(^)(0 = ll[uf\l)/wf\l)] = r(o)^-eFi.)d/ksT ^ ^ (18) 
J=l 

Solving this determines xs — Isd and thence yields the measured spring or stalling 
force 

= kBT{lnT)/9d = kBT£/9d, (19) 

where we have, of course, identified the zero-load rate factor, r*^°\ with the original 
rate factor F in Eq. 2. 
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It is striking that this expression for the staUing force (which depends on the 
quasiequihbrium hypothesis, Eq. 16, that is needed to extend the original kinetic 
model) agrees precisely with the barometric expression Eq. 13 for fs, provided 
one accepts the natural, near-equilibrium evaluation 9 — 1. We regard this overall 
consistency as strengthening both approaches. 

Velocity versus Load. The extended rate constants Uj and Wj given in 
Eq. 16 also serve to provide a relation for the motor velocity, V{F), as a function 
of a steady load force, F [and, equally, for the load-dependent diffusion constant, 
D{F)]. For arbitrary N one may appeal to Eq. Al which shows, as expected, that 
the stalling load, Fs, which brings V{F) to zero, agrees with Eq. 19, i.e., Fs = fs- To 
write an explicit result for iV = 2 in an illuminating form, we introduce the reduced 
force and modified load factors 

V^F/Fs and Af^^-{ef/e). (20) 

Then from Eqs. 3, 16 and 19 we obtain 

V{F) ^ asinh[|s(l-r^)]/sinh(|£) 

^(0) me-^^^^'' + U2e-^i^'' + wie^2^v ^yj^e^^ev' 

where, naturally, V{0) is simply the no-load result of Eq. 3; thus the right hand side 
reduces to unity when rj = (and vanishes as — >• 1). 

Now for £ small (say, < 2), so that the motor is operating not too far from 
equihbrium, one has V{F) fa T^(0)(1 — rj)/ {1 + cerj). This represents a hyperbolic 
force law which will be concave or convex depending on the sign, -|- or — , of c : see 
the illustrative examples in Fig. 1. For small c the law is close to linear and, in 
fact, c vanishes whenever Ui/S.2 + U2At = '"^1^2 + W2Ai- This condition has many 
solutions; for example, if the backward rates are small, so that 5 = {W1+W2) / {U1+U2) 
< 0.1, say, the loading scheme ~ fs ^9/ (1 -|- S) yields a near-vanishing c. 

On the other hand, if Ui greatly exceeds U2, Wi, and W2, the reduced {V, F) plots 
become insensitive to Ui. Then if, as mentioned (6), one has Ui ~ A;i[ATP], the plots 
will become independent of the ATP concentration (6). Furthermore, if £ is large but 
{92 /9)£ ~ 1, the {V, F) plots will be close to linear. 

Although straight, convex, and concave velocity-load plots are readily generated, 
other reasonable values of the six parameters: e, Wi/ui, W2/U1, and 92/9, 9^/9 and 
92 19, yield plots exhibiting points of inflection of either sense, as shown in Fig. 1. 
However, plots with negative inflection points, such as (e), are realized in relatively 
small regions of the parameter space. If negative 92 or 9^ are admitted (see after Eq. 
17) the velocity may initially rise when a load is imposed! Plots with two inflection 
points arc then also allowed. Thus if one could determine plausible values for the 
no-load transition-rate ratios, experimental {V, F) plots might, at least for an N — 2 
model, throw some light on the load distribution parameters, 9f. 
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Figure 1: Examples of velocity-load plots for = 2 models with various parameter 

sets {e; (wi, W2)/ui] {9+, 9^ = 92)/9}: (a) {lO^^; 0.99, 0.99; i, 0}, 

(b) {9.2; 10-2, 10-2; 1,0}, (c) {10-^; 10-^, 10-^; 0, i}, (d) {11.1; 10-^ 0.15; 0, i}, 

(e) {23.0; 10-^ 10-^ 0.07, 0.43}. Note that ^(0) is the velocity at zero load (Eq. 3) 

while Fs denotes the stalling load. 



Discrete Jump Models. As mentioned above, the Einstein force scale ob- 
tained from the kinetic scheme, Eq. 1, is subject to a fairly stringent bound. Analyz- 
ing the expression fE/ksT = V/D (see Eq. 8) for the case N — 2 (using Eqs. 3-6) 
we can prove that for all rates Uj, Wj > 0, the Einstein scale satisfies fE/kBT < 2N/d; 
the maximum for A'^ = 2 is realized for uniform rates Uj = mq ^ Wj = wq (all j), and 
we believe the same condition yields the bound as stated generally for all A^. (The 
uniformity condition can be understood heuristically since in such a case there arc 
no distinguishing rate-limiting steps in the cycle. The N — 2 model studied in (17) 
also respects the /ou;er bound fE/kBT > 2/d; but this is attributable to the special 
limiting situation, W2 = 0, studied there which cannot be literally true.) 

Our purpose here is to demonstrate that these bounds on Je are related to the 
continuous-time picture of the rate process embodied in the kinetic master equations 
based on Eq. 1; in essence, these force a minimum value of the diffusion constant D. 
To see this most directly, consider an {N — l)-state model with master equation 

dP 

-^{l, t) = uPo{l -l,t)+ wPo{l + l,t)-{u + w)Po{l, t), (22) 
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where we have put ui = u > wi = w > 0. Then one finds 

V^{u-w)d, D^^{u + w)(f, (23) 

[see, e.g. (16)]. Note the lower limit D > ^u(P, which is approached when w/u ^ 0. 
This leads directly to the bound fE/ksT < 2/d. 

By contrast, consider a discrete event sequence in which a forward or backward 
jump is attempted at (mean) time intervals At = r (triggered, one might picture 
for a molecular motor, by the arrival of individual ATP molecules). If Po{l;n) is 
the probability that the (motor) particle is at site / after n jump attempts, one has 
(19,20,24) 

Po{l; n + 1) = p+Po{l - 1; n) + poPo{l; n) + p_Po{l + 1; n), (24) 

where p+ and p_ are the probabilities of completing a positive or negative step while 
Po = 1 — P+ — P- is the probability of remaining at the same site. If one sets p+ = ut 
and p_ = WT, and identifies the time as t ~ nr, this discrete master equation, reduces 
to the continuous form, Eq. 22, in the limit r (24). 

Now the mean displacement {x)n after n = 1 attempts is clearly {p^ — p-)d so 
that the mean velocity is 

V ^{p+- P-)d/T ={u- w)d. (25) 

Note that the identifications appropriate to the continuous limit yield agreement with 
Eq. 23. To compute D we may use Eq. 7 with only a short time, interval, specifically 
t = T, since, by assumption, successive jump attempts are uncorrelated. Thus, from 

(a;^)i = (p+ + p-)d'^ we obtain 

D = i(dVr)[p++p_-(p+-p_)^] 

= l[u + w-(u-w)^T]d^. (26) 

To see that D now has no positive lower bound, we may specialize to the case 
Po = or consider the limit p-{— wt) <^ p+: then one finds D oc (rf^/r)p_|_(l — p+) 
which becomes indefinitely small when p+ approaches unity. Hence no upper bound 
on /e exists in such a discrete jump model. Indeed, it is intuitively clear that in the 
limit p+ = 1 (so that pq = p_ = 0) the particle moves essentially ballistically at speed 
d/r with no dispersion. 

Notice also that the barometric formulation can be applied directly to the jump 
model by using Eq. 24. It leads precisely to the previous form, Eq. 13, but with 
r = however, this agrees exactly with the continuous-time (iV = 1) result 

r = u/w when, as above, one puts p+ = ut and p_ = wt. In addition, the ratio 
R{s) = Jb/ Je obeys Eq. 14 but with, in leading order, (w/cr) replaced by = 
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(1 — po)/2(r + 1). For £ > 2 one has < 0.06 and R{e) varies much as discussed 
above for the continuous case. 

It might be objected that our arguments have more or less assumed that the 
jump attempts occur regularly at times nr whereas, more realistically, there should 
be some distribution, say ip{t), of waiting times between one event and the next. 
Then r would be the mean time between attempts, defined by 



Such a model may be studied along the lines of Montroll and Scher (24). Provided 
ilj{t) decreases sufficiently fast when t — > oo that the second moment is finite, the 
analysis for V and D can be carried through: it shows again that D is unbounded 
below while Je oc V/D is unbounded above. Indeed, Eq. 25 for V remains valid 
while Eq. 26 for D gains a factor (1 — 0) before each squared term, {p^ — pS)'^ and 
{u — wY. The parameter 9 = (t'^ — t^)/f^ > measures the relative width or "spread" 
of the waiting time distribution ipit)'- e.g., for ip{t) oc t^~^e~'^* with i/, 7 > 0, one hasQ 
r = v/'-f and O = l/p. The sharp distribution originally envisaged corresponds to 
the limit u 00. 

Finally, note that we can also analyze precisely multistate versions of these discrete 
jump models with waiting times. 



In order to understand the driving force, /, exerted by a molecular motor that takes 
steps of size d on a molecular track, we have analyzed a broad class of stochastic 
models: in particular, Eqs. 1 and 2 embody a general, "linear" reaction sequence. In 
the presence of a constant free energy source, the motor will achieve a steady velocity 
V (> 0) but with fiuctuations about the mean position, {x{t)) = Vt, described by a 
diffusion constant, D, that may be measured by observing the variance: see Eq. 7. 
The table below lists various force scales that then arise and their relation to /. 

By way of a concrete numerical illustration, consider a kinesin molecule moving on 
a microtubule (2-4, 6-10, 12) for which d ~ 8.2 nm (9). Svoboda, Mitra, and Block 
(6) observed V ~ 670 nm/s when [ATP] = 2 mM and measured the variance from 
which we obtain D ~ 1395 nmVs. At T = 300 K these data yield Je ^ 2.0 pN. On 
the other hand, the observed stalling force was fs — 5-6 pN (6, 7) significantly larger 
than Je, as we have argued it should be. Note also, comparing with the maximal 
force, fmax — 6.2 pN, that the observed efficiency £, is in the range 80-95%. (This 

^Thc specific results quoted in Rcf. (24), Eqs. 75 for = ^ and ly — 2 are in error (and the 
factor 4 in Eq. 76 sliould read 2). 



T = t with 




(27) 
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Table 1: Forces related to a MoleculEir Motor 



Maximum driving force 


fmax 


= AGo/d > f 


[Eq. 6] 


Einstein Scale 


fE-- 


= kBTV/D < f 


[Eq. 8] 


Gravitational force 


fa- 


= mg 


[Eq. 12] 


Barometric force 


fB- 


= ksTK = f 


[Eqs. 12,13] 


Stalling force 


fs = 


-- kBTe/Od = f 


[Eq. 19] 


Load and stalling load 


F,Fs, 


V = F/Fs 


[Eq. 20] 



observational estimate does not allow for the possible "wastage" of ATP by futile 
hydrolysis (13) without translocation of the motor: recall the footnote below Eq. 5.) 

The "barometric" force scale, /b, arises by considering an obstacle that blocks the 
motor's motion on the track: the resulting statistically stationary distribution decays 
with the distance z from the obstacle as e~'^^: see Eq. 9. It would be interesting 
(although difficult) to measure k, and to compare fs, so derived, with the observed 
stalling force fs- 

For the general {N — 2)-state model (see Eqs. 1-5) with transition rates ui, U2, 
Wi and W2 one has = (ksT/d) x \n{uiU2/wiW2)- For kinesin (from Drosophila) 
Gilbert and Johnson (12) studied the kinetics using chemical-quench flow methods. 
Assuming [ATP] = 2 mM their data show that Ui = 3800 s~^, U2 = 15 s^\ and 
Wi — 200 s~^ represents a sensible map on to an = 2 model; however, W2 proved 
unobservably small. Merely for illustration, therefore, suppose W2 — 1*2/100 = 0.15 
s~^. This gives ~ 116 nm/s and D ~ 474 nm^/s, which yield /e — 1-0 pN (at 
T = 300 K), while the rates give Jb — 3.8 pN. The agreement with the results of 
Svoboda et al. is not impressive: nevertheless, the orders of magnitude, the inequality 
Ib > Je, and the rough equality Jb — fs, are in full accord with our analysis. 

More recently, Higuchi et al (9) obtained data (for bovine brain kinesin) leading 
us to Ui ~ 1400 s~^ and U2 — 45 s"-*^, in only rough agreement with the values 
derived from (12). The ad hoc assumption wi/ui ~ W2/U2 — 1/100 yields V ~ 354 
nm/s and D ~ 1370 nm^/s, closer to observations (6). Likewise, Je — 1-1 pN and 
/b ~ 4.7 pN, now accord better with the direct experiments (although depending 
logarithmically on wi and ^2). While the general theoretical picture is supported, 
further experiments on standardized kinesin samples would clearly be valuable and 
could provide more stringent tests. 

To discuss the velocity V(yF) of a motor under a load F, the transition rates must 
be modified: see Eqs. 16 and 17 where the load-distribution factors, df, recognize 
that the various transitions in a motor protein probably accept quite different frac- 
tions of the total stress. Indeed, some forward rates might even be accelerated which 
could provide a mechanism to conserve, e.g., ATP under "no-load" conditions. It is 
natural to take the overall load-distribution factor 9 (Eq. 17) as unity which leads 
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to the equality of fs and fs- see Eq. 19 et seq. However, 6=1 can be doubted for 
real motors and might well be tested by experiment or simulation. 

Even for a two-state model, the expression for V{F) is quite complex: see Eq. 
21. As seen in Fig. 1, the six independent parameters permit velocity-load plots of 
varied shapes (including nonmonotonic forms not shown). Certain types, such as (e) 
characterize small regions of the parameter space; but, in general, the variation of 
V with F may reveal comparatively little about the motor mechanism or parameter 
values. 

Negative, i.e., assisting loads {F < 0) are predicted to speed up the motor and this 
has been observed (10). Conversely, under super-stalling loads {F > Fg), backwards 
velocities are predicted; single reverse steps of kinesin have then been seen (10) but 
no steady reverse velocities have been reported. This probably reflects very small 
terminal reverse rates, (12). Indeed, these transitions presumably describe second 
(or higher) order chemical reactions controlled by the low concentrations of hydrolysis 
products. The frequently observed process of detachment from the track (6,10) should 
also be included in a fuller account. 

The adequacy of the stochastic models encompassed by Eq. 1 is challenged by 
lower bounds on the diffusion constant, D, which yield upper bounds on /g. For 
kinesin at T = 300 K this bound is 2.03 pN for any (N = 2)-state model. The data 
of Svoboda et al. (6,7) essentially meet this; but were the bound violated, one might 
conclude that an A*" > 3 kinetic model was needed. However, models in which the 
transitions are described by discrete jumps occurring after certain waiting times, are 
not susceptible to these constraints. Such models might well prove more realistic, 
although at present the simpler kinetic representations seem adequate. Nevertheless, 
it should be noted that the main principles we have enunciated are not restricted 
to the N = 2 sequential kinetic models specifically analyzed. Consequently, the 
observation of significant violations would indicate serious deficiencies in the general 
understanding of motor mechanisms. 



We are indebted to David A. Huse, Stanislas Leibler, Michelle D. Wang and Benjamin 
Widom for valuable comments on our work. The support of the National Science 
Foundation (under Grant CHE 96-14495) is gratefully acknowledged. 

Appendix 

For a one-dimensional hopping model with N states and arbitrary transition rates 
Uj and Wj, as introduced in Eq. 1, Derrida (16) obtained the exact steady-state 
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behavior. For the drift velocity he found 

d 



N-l 



V 



j=0 % 



(Al) 



where d is the spatial period (or step-size) while 

1 / N-l k 



w 



k=l i=l "'j+i 



(A2) 



The expression for the diffusion constant (16) is similar but more complex and less 
illuminating. 

For N = 2, solutions can be obtained for all times following the procedure outlined 
in (17). Thus, the probability that the particle is at site I in state ji {j — 0, 1) after 
a time t having started at the origin, / = 0, is 



-iq{l+j/2) 



E+{q)e 



A+(?)t _ ^ 



5-(?)e 



A_(?)t 



with rate parameters X+{q) and X-{q) given by 



A±(g) 



-a ± ^Ja'^ + 4uiU2{e'^'i - 1) + 4wiW2(e-2i'? - 1) 



(A3) 



(A4) 



while a is in Eq. 5 and the coefficient functions are 



^^^^^ - X-{Q)-K{q) [ ^26^' + ^le-'' ) ■ ^ ' 
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